Fast growing instabilities for non-parallel flows 
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Unstable modes growing when two plasma shells cross over a background plasma at arbitrary 
angle 8, are investigated using a non-relativistic three cold fluids model. Parallel flows with 8 — 
are slightly more unstable than anti-parallel ones with 8 — n. The case 8 = tv/2 is as unstable as 
the 8 = one, but the fastest growing modes are oblique. While the most unstable wave vector vary 
with orientation, its growth rate slightly evolves and there is no such thing as a stable configuration. 
A number of exact results can be derived, especially for the 8 = n/2 case. 

PACS numbers: 52.35.Qz - 52.40.Mj 



I. INTRODUCTION 

Beam plasma instabilities are ubiquitous in physics and have been investigated for many decades The topic 
currently undergoes a renewed interest through the Fast Ignition Scenario for Inertial Fusion [2j or some scenarios of 
Gamma Ray Bursts production in astrophysics 3] . Indeed, Astrophysics offers a very wide range of unstable systems 
which can be magnetized, relativistic, homogenous or not. Counter-streams instabilities are usually studied assuming 
parallel streams. But why should real systems systematically fit this scheme? Admittedly, when an electron beam 
enters a plasma, a return current is prompted in opposite direction to neutralize it [![. But two already current and 
charge neutralized plasma shells could perfectly collide over a background plasma at an arbitrary angle. It seems the 
problem of non-parallel streams has not be addressed so far, and the goal of this letter is to show that non negligible 
growth rates can arise from non-parallel streams interactions as well. 

Let us consider two non-relativistic beams, both charge and current neutralized. For simplicity, we consider here two 
electron-proton beams where both species have equal densities so that charge and current neutrality are guarantied 
regardless of the beams relative motion. The problem of the beams respective orientation is interesting only if there 
is a background plasma. Otherwise, one just needs to consider the reference frame of one of the beams to cancel any 
orientation parameter. We eventually come up with the setting pictured of Figure [TJ Calculations are conducted in 
the reference frame where the background plasma of electronic and protonic density N is at rest. In order to focus 
on the parameter 9, the two beams have equal protonic and electronic densities n, and both flow at the same velocity 
V along their respective direction. We now study the stability of harmonic perturbations oc exp(ik • r — iujt) where 



X Background, N p = N e = N 



Beam 2, V 2 =V 
n p 2 = n p 2 = n 



j 9 Beaml,Vi=V 
^ n p l = n p i = n 



FIG. 1: Scenario considered. 
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FIG. 2: (Color Online) Growth rate map in terms of the reduced wave vector Z for a = 0.1, f3 = 0.1 and 8 — 0, n/2 and n. 



i 2 = — 1 and k = (k x , 0, fc z ). Note that the all k unstable spectrum is evaluated in order to be able to spot the fastest 
growing modes for any given configuration. Finally, we neglect proton velocities perturbations in view of their much 
larger inertia. 



II. THREE FLUIDS MODEL 



As a first approximation, we implement a cold three fluids model for the two beams and the plasma. Linearizing 
the conservation and Euler equations for the three species, the dispersion equation is found evaluating the dispersion 
tensor in a standard way in terms of the dimcnsionlcss variables, 

a =T7! Z = P = —, X= , (1) 

ly ojp c ujp 

where to 2 — AirNe 2 /m e is the background plasma frequency, e the electron charge and m the electron mass. Due 
to the arbitrary orientation of the wave vector and of the two beams, the full dispersion tensor is too large to be 
reported here. Still, the dispersion equation remains polynomial and can easily be solved numerically. For Z x = 0, 
the zz component of the tensor give the growth rate of a two-stream like instability, as least for 9 = and n. The 
corresponding dispersion equation reads, 

la a 

1_ ^2 - ( X -Z z ) 2 ~ {x-Z z cos9) 2 = ' ( ' 

The growth rate map in terms of Z for 9 = 0, ir/2 and it on Figured] The first and the last cases pertain to 
well-known systems as 9 = eventually comes down to one single beam of density 2n interacting with the plasma, 
while 9 — 7T corresponds to two counter-streams crossing over a background plasma. These situations have been well 
studied, and it is known that within the present non-relativistic regime [5(, they are governed by the two-stream 
instability which dispersion equation is precisely given by Eq. @. The maximum growth rate 8m is in the diluted 
beam regime a«l, 

6 M (d = 0) ~ ^a 1 / 3 ^, and 5 M (6 = tt) ~ ^a 1 ^. (3) 

For 9 = 0, the two beams act as one, and a resonant unstable mode is feed by the free energy of both at the same 
time. For 9 = n, resonant unstable modes can only travel with one single beam, yielding a smaller growth rate. 



III. MAXIMUM GROWTH RATE FOR 9 = vr/2 

Besides these extreme orientations of the two beams, Fig. [5] clearly displays some interesting features for the case 
9 = tt/2. One the one hand, the fastest growing mode is here found for Z z = Z x ~ 1 so that it can stay is phase and 
exchange energy with both beams at the same time. Such non-trivial orientation of the fast growing mode has been so 
far related to relativistic effects [H, 0] ■ We find here that some unusual system geometry can produce the same effect as 
the only way for a mode to move in phase with two non-parallel beams is to follow an oblique direction. On the other 



3 



5 

0.3 
0.2 
0.1 



V 



3jr 
4 



a = 0.1 

-a = 0.01 
a = 0.001 

e 



FIG. 3: Evolution of the maximum growth rate (lu p units) for j3 == 0.1, in terms of the angle 9 and for various beam to plasma 
density ratios a. 



hand, the maximum growth rate seems very close to the one reached for 9 = 0, and some finer numerical evaluation 
shows that it is almost the same. Such equality can be demonstrated setting 9 = n/2 from the very beginning, and 
considering the growth rate for wave vectors fulfilling Z z = Z x . The dispersion equation for these modes reads, 



ar [x 



2a 



P{x) = 0, 



(4) 



where 



P(x) = {x 2 - l)(x - Z z f - 2a(Z 2 z +(x- Z z f) - 2(x - Z z ) 2 Z 2 /f3 2 . (5) 

Since the first factor in Eq. (H} is strictly equal to the dispersion equation @ for 9 = 0, we here prove that the 
maximum growth rate for this branch is rigourously the same, and that it is reached for the very same Z z component. 
A closer look at P{x) shows that it yields another unstable mode. By setting Z z ~ 1 in the expression of P(x) and 
developing the result near x = 1, we find the growth rate of this second oblique unstable mode, 

hi2 ~ uJpPVa- (6) 
For a and (3 lower than unity, this secondary growth rate remains smaller than 8 m (9 = 0). 



IV. CONCLUSION 



The overall system is thus found exactly as unstable for 9 — as it is for 9 = ir/2. Beyond this value of 9, the 
maximum growth rate must decrease to its final value 8m (6 = n), as given by Eq. ([3]). Figure [3] shows how 8m evolves 
between 9 = and n, for (3 = 0.1 and various beam to plasma density ratios a. As expected, 8m remains almost 
constant between and 7r/2 while it falls down to its final value <5m(t) for 9 ~ 2. Because 8m is independent of (3 for 
9 = 0, 7r/2 and it, we can expect an overall weak dependance for other angles. This is confirmed through numerical 
calculation as the three curved displayed on Fig. [3] are indistinguishable from their (3 = 0.01 counterparts. 

Ton conclude, it is important to emphasize that there is no such thing as a stable configuration. On the contrary, 
the system remains unstable regardless of the beams orientation, and the evolution of the maximum growth rate is 
limited. The most unstable wave vector is two-stream like for parallel and anti-parallel orientations, but turns oblique 
in the intermediate case to stay in phase with both beams at the same time. Noteworthy, relativistic and/or kinetic 
effects should force an oblique regime regardless of the beam orientation so that an evaluation of the whole unstable 
spectrum becomes mandatory. 
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